Mechanics and force transmission in soft composites of rods in elastic gels 
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We report detailed theoretical investigations of the micro- mechanics and bulk elastic properties of 
composites consisting of randomly distributed stiff fibers embedded in an elastic matrix in two and 
three dimensions. Recent experiments published in Physical Review Letters [102, 188303 (2009)] 
have suggested that the inclusion of stiff microtubules in a softer, nearly incompressible biopolymer 
matrix can lead to emergent compressibility. This can be understood in terms of the enhancement of 
the compressibility of the composite relative to its shear compliance as a result of the addition of stiff 
rod-hke inclusions. We show that the Poisson's ratio u of such a composite evolves with increasing 
rod density towards a particular value, or fixed point, independent of the material properties of 
the matrix, so long as it has a finite initial compressibility. This fixed point is v = 1/4 in three 
dimensions and = 1/3 in two dimensions. Our results suggest an important role for stiff filaments 
such as microtubules and stress fibers in cell mechanics. At the same time, our work has a wider 
elasticity context, with potential applications to composite elastic media with a wide separation of 
scales in stiffness of its constituents such as carbon nanotube-polymer composites, which have been 
shown to have highly tunable mechanics. 

PACS numbers: 87.16.Ka, 62.20.Dc, 82.35.Pq 



I. INTRODUCTION 

Many materials, natural (T], 0] as well as man-made 
Q, that we come across in our daily lives are composite 
materials, combining multiple components with distinct 
elastic properties. Notable examples of natural compos- 
ites are wood, bone, plant and animal cells, while glass 
or carbon fibers embedded in epoxy resins are two exten- 
sively used man-made composites. A hallmark of such 
materials is that the components they are made up of 
interact in a highly synergistic manner such that the col- 
lective properties are more than merely the sum total 
of those of the constituents. These composites are often 
made up of a soft elastic background, reinforced with stiff 
fibers. Because of their aspect ratio, fiber reinforcements 
have inherent advantages over other geometries allowing, 
for example, for long range force propagation or enhanced 
strength along a given direction. By tuning the concen- 
tration of the fibers and the relative mechanical proper- 
ties of the fibers and the matrix, one can modulate the 
bulk properties of the composite resulting in materials 
with remarkable properties. Fiber-reinforced composites 
often form the building block of structures that require 
high strength and stiffness as well as low weight, both 
in materials optimized by nature through evolutionary 
processes as well as those engineered by humans. 

A prime example of a natural or living composite is the 
cell cytoskeleton, a composite polymeric scaffold made up 
of several distinct filamentous proteins of varying lengths 
and stiffnesses, which helps the cell maintain its shape 
and provides it support. Most previous biophysical stud- 
ies of cytoskeletal networks have focused on purified gels 
or networks consisting of one type of filament [3-l22| . The 
cytoskeleton, however, contains three major types of fila- 



ments: microtubules (MTs), filamentous actin (F-actin), 
and intermediate filaments. F-actin and intermediate fil- 
aments show thermal bending on micron length scales, 
and often act collectively as networks with a sub-micron 
mesh size. MTs, on the other hand, have a persistence 
length of the order of millimeters and effectively behave 
as rigid rods on cellular lengthscales. The mechanics of 
a cytoskeletal composite consisting of F-actin and MTs 
can, therefore, be better understood by modeling it as 
a fiber-reinforced elastic composite with stiff rods ran- 
domly and isotropically distributed in a much softer elas- 
tic matrix, rather than in terms of a network or elastic 
continuum made of just one type of material. There have 
been many fundamental studies of the engineering prop- 
erties of composite materials [23l - l28j that attest to the 
highly synergistic collective material properties of these 
materials. However, these studies often assume compo- 
nents of comparable stiffnesses or fibers aligned in special 
directions, and hence have limited applicability to the bi- 
ological composites discussed above. 

Recent studies [29'-'33!| on reconstituted composite cy- 
toskeletal networks have shown that synergistic mechan- 
ical interaction between F-actin and microtubules indeed 
lead to viscoelastic properties very distinct from one- 
component networks |4l422l| . A remarkable recent experi- 
mental finding on the elasticity of composite cytoskeletal 
networks was the appearance of enhanced compressibil- 
ity of F-actin networks upon the addition of MTs p^. 
Many materials found in nature are nearly incompress- 
ible, which implies that when under external mechanical 
stress or strain, they try to conserve their volume, by 
changing their shape. This is because they have a much 
larger bulk modulus than shear modulus. It is intrigu- 
ing, therefore, that the addition of very stiff rod-like MTs 
would give rise to a finite compressibility for such a com- 
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posite of F-actin and microtubules, while pure F-actin 
matrices appear to be incompressible. 

In this paper we report detailed investigations on the 
mechanical response of a composite of MTs embedded in 
an F-actin matrix by modeling it as a composite mate- 
rial consisting of rods in an elastic matrix [3J] . We use a 
dipole approximation for the rod-like inclusions and treat 
the composite as an effective medium whose properties 
we study using both micro-mechanical and continuum 
approaches. We carry out our investigations for elastic 
matrices in both three and two dimensions. We report 
flow diagrams for the Poisson's ratio v with increasing 
rod density that exhibit a stable fixed point at = 1/4 
in 3D and at = 1/3 in 2D, both of which correspond 
to a Cauchy solid with equal Lame coefficients A and 
Thus, adding rods to 3D matrices characterized by a 
Poisson's ratio 1/4 < < 1/2 makes the material more 
compressible relative to the shear compliance (i.e., u de- 
creases), while for materials with v < 1/4, stiff rods lead 
to a less compressible medium. In 2D, we find enhanced 
compressibility for matrices characterized by Poisson's 
ratio 1/3 < < 1, and a decrease in compressibility for 
those with a Poisson's ratio v < 1/3. We also evaluate 
the effective medium elastic moduli of the composite as 
functions of the concentration of rod-like inclusions using 
a self-consistent approximation. While our self-consistent 
approach is only approximate at intermediate concentra- 
tions, we obtain exact results in the limits of high and 
low concentrations. 

The rest of this paper is organized as follows. We de- 
scribe our model and methods of exploring the collective 
mechanics of the composite under investigation in Section 
im in III Al we describe a continuum approach to study- 
ing a fiber-reinforced composite with a 3D matrix, while 
in section III Bl we elucidate a micromechanical method 
informed by linear response theory. We discuss the self- 
consistent calculation and the ensuing results in Section 
mil In section IIVI we describe an approach for study- 
ing this system that takes into account the full tension 
profile along the fibers. In section |V] we describe the 
calculations for a two dimensional matrix. We conclude 
with a discussion of the implications of our findings for 
the cytoskeleton and other composite materials in Sec- 
tion ED 



II. MODEL 

The main focus of our study is a composite of F-actin 
and MTs in the cell cytoskeleton. As discussed earlier, 
MTs are three orders of magnitude stiffer than F-actin. 
F-actin forms a three dimensional meshwork in the cy- 
toskeleton, with an average mesh size of ~ lOOnm. MTs, 
on the other hand, effectively behave as rigid rods and 
are generally 1 — 10 microns long in most cells, whereas 
in axons their length can be 50 — 100 microns. The 
simplest model that can capture the essential underly- 
ing mechanics of this composite is then a collection of 




FIG. 1: (Color Online) We consider a representative element 
of the effective medium made of the background elastic matrix 
and rods. Here we show one representative rod oriented in a 
direction as shown, with a polar angle 6 and azimuthal angle 



isotropically distributed stiff elastic rods in a compara- 
tively soft elastic background. We further consider the 
limit where there is no sliding of the MTs relative to the 
matrix. There is increasing experimental evidence that 
MTs are directly mechanically coupled to other struc- 
tures in the cytoplasm, including F-actin matrices. The 
nature of this interaction is not fully understood and it 
has been proposed that these interactions are mediated 
by microtubule associating proteins (MAPS). In the limit 
that the strength of this coupling is large, and the move- 
ment of the MTs is further stearically hindered by the 
presence of other stiff structural elements, such as other 
MTs or stress fibers, the MTs can be considered rigidly 
embedded in their matrix. We construct an isotropic and 
homogeneous effective medium with the same mechanics 
as this fiber reinforced composite of stiff rods embed- 
ded in an elastic background. We calculate the effective 
medium shear modulus /i and longitudinal modulus A 
for our system, followed by other mechanical quantities 
as such the Poisson's ratio. We use two different methods 
in our study: a macroscopic approach guided by methods 
in continuum mechanics, and a micro-mechanical method 
based on linear response theory. 



A. Continuum Approach 

For an isotropic and homogeneous elastic material un- 
der the action of external forces, the stress tensor aij 
is related to the strain tensor Uij by the expression 
(Ty = XSijUii + 2^Uij, where fj, and A are the Lame con- 
stants of the material, 5ij is the Kronecker Delta func- 
tion and summation over repeated indices is implied. 
For such a material, the displacement field Ui at a po- 
sition r due to a force / acting at point r' is given by 
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Ui{r) = aij{f—r')fj{r'), where is the elastic response 
function or Green's function for the material. The re- 
sponse function aij (r) has just two distinct components, 
corresponding to the response parallel and perpendicular 
to r: 
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[r,f,(l-/?)+%(l + /?)], 



(1) 



where 



1 



(2) 



Here, /3 provides a simple measure of the degree of com- 
pressibility of the material, relative to the shear compli- 
ance. Specifically, it is given by the ratio of the shear 
modulus /Lt to the longitudinal modulus: 



A + 2^ 



(3) 



For incompressible materials, (3 — and the parallel and 
perpendicular response functions are related by a simple 
factor of two: a^^{r) = 2aj_{r) = l/An^r. These corre- 
spond to the elastic analogue of the Oseen Tensor [s^ 
for incompressible materials in 3D. 

A related measure of a material's compressibility is 
the Poisson's ratio. It is the ratio, when a material is 
stretched or compressed in one direction due to an ex- 
ternally applied force, of the transverse strain (perpen- 
dicular to the applied force), to the axial strain (in the 
direction of the applied force) . The Poisson's ratio is de- 
fined in terms of the Lame constants as u = ^A/ (A -I- /i) 
in 3D and has an upper bound of 0.5 (for incompressible 
materials) and a lower bound of —1 (corresponding to 
zero bulk modulus). For a composite material made up 
of an elastic matrix and rods, the addition of rods will 
lead to additional stresses in the medium and changes in 
the Lame coefficients, compressibility and Poisson's ratio 
of the composite, which we calculate as follows. 

We model the composite as an isotropic and homoge- 
neous effective medium that is made of the bare elas- 
tic medium (e.g., the F-actin matrix) and a collection of 
rods (MTs) embedded in it, and with elastic properties 
macroscopically indistinguishable from those of the ac- 
tual composite. We consider a 3D sample of this material 
under uniaxial compression or extension, i.e., we apply a 
uniaxial strain tzz along the z direction, with all other 
strain components zero. We consider a rod of length a 
embedded in this effective medium, making polar and 
azimuthal angles 9 and (j) with respect to the coordinate 
axes as shown in the schematic Fig. [T] When the ma- 
terial is strained as described above, the rod suffers a 
change ea in end-to-end distance, where e = tzz cos^ 9. 
For the sake of simplicity, we use a dipole approximation 
for the constraint of fixed length a of the rods that is 



assumed to be small compared to all other length scales 
in the problem. The strength of the induced dipole can 
be calculated as follows. Replacing the rod by a (tensile) 
point dipole of strength p and orientation a at the center 
of mass of the rod, we obtain a net change in end-to-end 
distance, which we set to zero. This gives 



= ea + 2pa\{a/2) 



ea 



2p_ 



(4) 



where the prime denotes derivative. The rigid rod is thus 
mechanically equivalent to a dipole of strength p = xe, 
where 



(5) 



The above represents the linear response of the medium. 
We note that there is also a linear order change in the 
orientation of the rod. Given that the strength of the 
induced dipole is linear in the strain, we can neglect the 
effect of the reorientation of the rod. 

For rod orientations in a given solid angle d£L — 
sin 9d9d(l), the stress arising due to these induced dipoles 
is given by (5 cry = ^xeaifij-dr^. The presence of these 
rods, with a number density n, will lead to changes 
in the Lame constants of the material, thereby giving 
rise to additional stresses at the boundary given by 
S(Txx = Sayy = SXezz and dazz — {26 + SX)ezz, for 
an isotropic distribution of such rods. Thus, the change 
in the Lame constants are given by: 

SX = nx f cos(6l)2sin(6')2cos(0)2 — 

J 477 



2S^l + SX = nx J cos{9) 
Solving for 6fi and SX, we find 



5^ ~ SX 
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Xn. 



(6) 



(7) 



Using the definition of /3 and Xj tti6 above expression 
translates to the equation 



d(3 = -a3/3(l - 3/3)n 



(8) 



for /?, the measure of compressibility relative to the shear 
compliance. Thus, /3 increases for < /3 < 1/3 and de- 
creases for (3 > 1/3, while it remains unchanged for (3 = 
and (3 = 1/3. This implies that for a medium that is only 
slightly compressible to begin with {/3 small, but posi- 
tive), adding rods makes it more compressible, while for 
a highly compressible medium {(3 > 1/3) the rods make 
it less so, as illustrated in Fig. [2l This suggests a stable 
fixed point (to the addition of rods) at /3 = 1/3 and hence 
Poisson's ratio v = 1/4. Similarly, (3 = and f = 1/2 
corresponds to an unstable fixed point. Our results are 
qualitatively consistent with recent microrheology exper- 
iments on a composite of microtubules embedded in fil- 
amentous actin 29] , which reported enhanced compress- 
ibility (v < 0.5) when stiff microtubules were added to an 
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FIG. 2: (Color Online) The flow diagram for the degree of 
compressibility /3 showing a stable flxed point at /3 = 1/3, 
and an unstable fixed point at = 0. 



almost incompressible actin matrix (i/ ~ 0.5), as inferred 
from the measured parallel and perpendicular response 
functions. 

So far we have considered the rods to be inextensible. 
We now consider the rods to have finite stretch modulus 
K = AEr, where A is the cross-sectional area of the 
rod with Young's modulus Er- We apply an extensional 
strain e and once again use a dipole approximation for the 
rod such that this dipole of strength p = can cancel 
the end-to-end displacement of the rod caused due to the 
applied strain. We can write down the resulting force 
balance for an extension A of the rod. 



givmg, 



X = 



_2p_ 



P_ 
K 



{2 + TTa^fi/K)' 



(9) 



(10) 



In our effective medium theory, the Lame coefficients 
A and /i for a composite with a density n of rods are, 
therefore, given by the following nonlinear relations: 



SX = Sfi = 



30 



(11) 



In the limit of very stiff rods, this reduces to the expres- 
sion derived earlier for inextensible rods, while for highly 
compliant rods, this is consistent with the elastic modu- 
lus Sfi = j^tfEj. of an affinely deforming rod network of 
volume fraction = Aan. 



B. Micromechanical method 

We now investigate the micro-mechanics of our system 
using linear response theory that also provides the guid- 
ing principles for the micro-rheology experiments used 
to study the F-actin-MT composite in (2^]. Here, we can 




FIG. 3: (Color Online) We consider a point-force fx applied 
at the origin (^4) and calculate the response at points B lo- 
cated parallel (here, along the a;— axis) and perpendicular 
(along the j/— axis) to the applied force. The rod center of 
mass is located at a distance r from the origin, and makes a 
polar angle 6i and azimuthal angle ifii. The rod is oriented in 
a direction as shown. 



think of this change in the response as arising from a 
cloud of induced dipoles in the elastic continuum. We 
calculate the displacement field Ui at a position r in an 
isotropic and homogeneous elastic material due to a force 
/ acting at the origin (for simplicity) , using the response 
function in Eq. ([T]). 

For a composite material made up of an elastic matrix 
and rods, the application of a force / at, say, a point A 
will elicit a net axial deformation of a rod located else- 
where in the medium due to its finite length. The pres- 
ence of rigid rods gives rise to constraints on the displace- 
ment field induced by the applied force. The collective 
elastic response of the composite at a position B depends 
on not only the applied force and the bare elasticity of 
the matrix, but also on the stiffness and concentration 
of rods. We calculate the change in the response func- 
tion and Lame coefficients upon addition of rods to the 
background elastic matrix by, once again, using a dipole 
approximation for the constraint of fixed length of the 
rods. 

We consider a single rod of length a and orienta- 
tion d embedded in the elastic medium as shown in the 
schematic figure [3l For simplicity, we consider an ap- 
plied force at the origin of magnitude / directed along 
the X-axis. The cm. of the rod is assumed to be at 
r = (x, y, z). We assume the rod length a is small com- 
pared with r and the separation between A and B. The 
force / leads to a net relative displacement of the ends 
of the rod, given by 



Au^ = (a ■ V)aij{r)Jj = lijkirjfjCik, 



(12) 
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and an axial strain, given by where p = r /h and 

e = ai{a ■ V)aij{f)fj = aaijkir)fjak, (13) 



where 



lijkir) = Vfca,j(f). 



(14) 



Here, we have kept only the leading terms in a. Also, 
we approximate the constraint of the rod by an induced 
dipole at its cm. of strength p = — X'^i'iklijk{r)fj- 
This gives rise to a displacement at a point h (represent- 
ing one of the points B in Fig. [3|): 

5ui(b) = pajCLkjijkib - r). (15) 

This defines a correction 

Soiij = Xlikiib ~ r)jj„in{r)akaid„idn (16) 

to the response function. For a uniform concentration 
n of isotropically distributed rods, the correction to the 
parallel component is given by 



foil (6) 



An 



Ixkiihx - f)jxmnir)dkdid,nd„ 



n / ((5aj|(&,p))p dp, 
Jo 



(17) 
I 



dndn' 

47r 



(18) 



xmn 

{bp)dkdidmd„ 



represents an average over the orientation O of a and 
an integral over the orientation il' of f. Similarly, by 
evaluating the displacement field at 6 = by, we obtain 
5a±ib). 



In order to evaluate Sa^\^j_, we perform both a full 
Taylor series expansion of (Sa\\ j_(b,p)) for small p (i.e., 
p < 1), as well as an asymptotic expansion for large p 
(i.e., p > 1). Following the angular integrals above, we 
are left with just three non-zero terms, in 0{p), 0(1/ p^) 
and 0(l/p^), giving the following simplified expressions 
for the change in the components of the response func- 
tion. 



p2((5a||(6,p)) 
p'{Sa^ib,p)) 



2{-3 + 213 + 11 13^) p 



450&V " 1 (-15 + 10/3 - 20/32)/p2 + (g _ 36/3 + 27 l3^)/p'^ 



Tia^ ^ j(3-2/3-11^2)p 



4506V ' ' [ (-15 + 10/3 - 20/32)/p2 + (ig _ 27/3 -I- 9/3^)/ p^ 

I 



if p< 1 
if p > 1, 

if p < 1 
if p > 1. 



(19) 



The remaining integrals over p lead to 
(5q;[| — — {tt/ 30) na^ a 
6aA_ = i(H-3/32)fo||. (20) 

Expressing a\\ and a± in terms of p and A using Eq. ^ 
and Eq. we obtain Sp = SX = j^X^ describing the 
change in the Lame coefficients of the effective medium 
with rod density, as in the continuum approach. 

III. SELF CONSISTENT CALCULATION 

We now use a self-consistent approximation wherein 
each added rod sees the composite as an isotropic and 
homogenous effective medium, with the Lame coeffi- 
cients A and p described by their effective medium val- 
ues. This is similar to self-consistent methods employed 



for aligned fiber-reinforced composites [23|. For a small 
increase in the number density of rods dn = n, the 
change in the compressibility /3, discussed in section 
II can now be described by the differential equation 
d(3/dn = ^a'^/3(l — 3(3). This suggests a way of calculat- 
ing the compressibility measure (3 and the Poisson's ratio 

V for composites with finite rod density. The resultant (3 
and V are shown in Fig. 3] for several different initial val- 
ues /3o in the absence of added rods. As earlier, we find 
that an incompressible material will stay incompressible 
(with Poisson's ratio v = 1/2) even on adding rods to 
it, but for a matrix with finite compressibility, however 
small, addition of rods tends to drive the system towards 
a state with (3 = 1/3 and hence v = 1/4. This suggests 
a stable fixed point (to the addition of rods) at /3 = 1/3 
and V = 1/4, and an unstable fixed point at /? = and 

V = 1/2. 

Further, this approximation allows us to cast Eq. (jlip 
as a set of differential equations representing the increase 
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FIG. 4: The degree of compressibility /3 and Poisson's ratio v 
of the composite as a function of mesh size ^ for inextensible 
rods, for different values of the bare degree of compressibility 
/3o of the medium. The solid (red) lines correspond to the 
stable and unstable fixed points at /3 = 1/3, = 1/4 and 
/? = 0, ^ = 1/2 respectively. The mesh size ^ is related to the 
rod density n by 1/^^ = na. 
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FIG. 5: (Color Online) The solid blue lines show the shear 
modulus of the composite as a function of the mesh size for 
different values of the ratio of the rod and medium compli- 
ance, while the dashed red line represents the affine result. 
The inset shows the Lame coefficient A for different values of 
the initial degree of compressibility /3o of the medium for both 
extensible and inextensible rods. The ratio ^o/ fir ~ 0.001 in 
the inset, except for the inextensible rod (dashed line), where 
/3o = 0.1. 



of the moduli upon the addition of stiff rods. While this 
represents an uncontrolled approximation and the result 
in Eq. pT|) is valid at small densities where the shear 
modulus fi on the right- hand- side can be approximated 
by that of the (bare) matrix, we find that integrating 
Eq. (jlip yields an exact expression in the limit of high 
density of the rods. The solution for fi is given by 



HrW \ — exp 



30 



(21) 



where, fiQ is the shear modulus of the medium in the 
absence of rods and = 2K/{na'^). Here, W{z) is the 
principal value of the Lambert W-function, which is de- 
fined by z = We^ . From this, we can also obtain the 
longitudinal modulus A as A = Aq + /i — /io once Aq and 
fiQ for the bare background matrix are known. The re- 
sults are shown in Fig. [5] for various initial conditions /xq 
and Aq. For small densities n and large K (such that 
Mo/^r ^ this reduces to ^ ~ fio{l + nna^/SO), 

since W{z) ~ z for small z. This is consistent with the 
above results for dilute systems with inextensible rods. 
As the density of rods and corresponding shear modulus 
increase, however, W{z) ~ ln(z). For highly compliant 
rods with fir ~ l^o, this is consistent with the elastic 
moduli of an affinely deforming rod network of volume 
fraction ip: Sfi = j^ipEr- Thus the addition of rods can 
significantly alter the collective elasticity of the compos- 
ite, and irrespective of whether the background elastic 
matrix is only slightly or considerably more compliant 
than the rods, we find that the composite stiffens signifi- 
cantly with added rods, finally approaching the elasticity 
of stiff affine rod networks at sufhcieirtly high rod densi- 
ties. 

Such a self-consistent approximation as employed here 



is only valid in the limit of strong direct interactions be- 
tween the rods. It is analogous to the self-consistent ho- 
mogenization approach introduced by Hill (23l | and fur- 
ther extended in differential effective medium theories 
[sj. The latter model two phase composites by incre- 
mentally adding inclusions of one phase until the desired 
proportion of constituents is reached. In the limit of 
weak/indirect interactions between the rods mediated by 
the surrounding matrix, such a self consistent approxima- 
tion may no longer be valid as each rod will then have 
a surrounding boundary layer (at length scales smaller 
than its length) where it probes the bare elasticity of the 
matrix and not of the composite. 



IV. CALCULATION FOR THE FULL TENSION 
PROFILE ALONG RODS 

In our calculations so far, we have used a dipole ap- 
proximation for the constraint on the rod. For an elastic 
rod of finite length, however, the displacement field varies 
more smoothly than for a dipole, decreasing as one goes 
from the end of the rod toward its center. Concomitantly 
the strain along the rod is uniform at its center and van- 
ishes at its ends, as shown in Fig. [Bl Our approximation 
of the extensional resistance of a rod by a dipole is thus 
expected to overestimate the contribution of the rod to 
the effective medium shear modulus. One can correctly 
account for the tension profile along the rod using an 
approach analogous to slender body theory in fluid dy- 
namics Q as follows. 

Let v{x) be the displacement field along the rod in the 
presence of a background strain e of the matrix. The ten- 
sion along the rod is given by Kv'{x), and the gradient of 
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FIG. 6: (Color Online) The contours of the tension Kv' (x) 
divided by e along a rod of length a = 1, as a function of 
distance x from the cm of the rod and the ratio \ir- 



this corresponds to a net force per unit length on the rod. 
This force is proportional to the relative displacement of 
that section of the rod with respect to the background 
medium. The displacement field can be obtained from 
the resulting condition for force balance: 



Kv"{x) = Cb(a;) 



(22) 



where C represents the elastic coupling of the rod to the 
matrix. It can be thought of as the drag coefficient per 
unit length of the rod, similar to the viscous drag on a 
slender body in the presence of a background velocity 
field [sHi- For an elastic medium, we approximate C = 
27r/i/ ln(^/c) where c is the rod cross-sectional radius, 
as in Ref. Here, the screening length ^ is of order 

the average separation or mesh size of the rod network, 
which varies with rod density. However, since C only has 
a weak logarithmic dependence on ^, we will treat it as a 
constant. Solving the above differential equation for v{x) 
using the condition of vanishing tension at the boundaries 
of the rod, Kv'(al2) = Kv'{-a/2) = 0, we have: 



v{x) 



£o sech 



24 



sinh 



(23) 



Here, £q — \J K/C, represents the length over which the 
longitudinal state of strain of the rod varies The 
tension profile along the rod is shown in Fig. [6] Once 
again, we can represent the rod as a force dipole on scales 
large compared with the rod length a, but with the dipole 
strength p = xe now given by p — ''"^^ ' 
therefore with x 



2Kv"{x)x dx, and 



X = Ka [1 - 2£o tanh (a/24) /a] 



(24) 



In the limit of highly compliant rods (small K)^ the 
lengthscale £q becomes a small length, and hence the 
strain and tension along the rod is nearly constant, ex- 
cept very close to the ends. For very stiff rods, on the 



FIG. 7: (Color Online) The dipole susceptibility x as a func- 
tion of K/ n for infinitely stiff rods (red dotted line), compliant 
rods with the dipole approximation (black dashed line) and 
compliant rods in the Batchelor-like calculation (blue solid 
line) with /x = 1, a = 1 and the factor In (^/c) set to 1/3. 



other hand, the strain exhibits a quadratic dependence, 
reaching a maximum at the center of the rod and vanish- 
ing at the ends of the rod. In this case, x = Cfl'^/12, which 
is, as expected, somewhat smaller than the value above 
for the simple dipole approximation for inextensible rods. 
Specifically, it is smaller by a factor of 31n(^/c). For in- 
termediate values of rod compliance, i.e., if one retains 
the 0(1/4^) term in the expansion of x given by Eq. 
([24| for small 1/4, we obtain x = Qa^il-Qa^ /1()K)/12. 
Thus, although the resulting differential equation for /i 
as a function of n is much more complicated if one takes 
into account the full tension profile along the rods, we 
note that, apart from the prefactor of 3 In (^/c) discussed 
above, the dipole strength derived in Eq. ^ takes on ex- 
actly the same limiting values for stiff and compliant rods 
as in Eq. (|24p . and approximates intermediate values to 
within no more than 13%. Figure[7]shows the comparison 
between x calculated above (Eq. (pi]) ) and for the case 
of inextensible and compliant rods calculated using the 
dipole approximation, as a function of rod stiffness. Fur- 
thermore, for rods that interact directly with each other, 
the elastic moduli of the composite can be obtained us- 
ing a self-consistent approximation [2^ as discussed in 
the previous section, while for rods that only interact 
with each other through their matrix, the rod concentra- 
tion only enters the calculation via the screening length 
^. Thus, the functional forms in Figs. |4] and [5] are 
expected to be good approximations. We find that both 
Lame coefficients once again evolve in the same way upon 
the addition of rods: dfi = dX = nx/15 This means 
that the qualitative form of dp/dn in Fig. [21 as well as 
our conclusions regarding the fixed points at /3 = and 
/? = 1/3, remain unchanged. 
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V. RODS IN A 2D ELASTIC MATRIX 

We now consider a planar random fiber composite, 
such as rods embedded in a two dimensional membrane. 
Examples of planar random fiber composites include or- 
dinary paper, synthetic and biological polymer mats and 



buckypaper [41[, the last one being a prime example of 
a planar composite that can bear large loads. We apply 
the methods described in the previous sections to study 
the elastic properties of such composites. Following the 
continuum method described earlier, we consider a 2D 
effective medium representing a planar rod-based com- 
posite in the x — y plane subject to uniaxial strain eyy 
along the y direction, with the constraint of zero strain 
along the x axis. We consider a single rod oriented at an 
angle 9 with respect to the y-axis. Let 6X and be the 
change in the Lame constants due to the addition of such 
rods, with density n and distributed isotropically, giving 
rise to extra stresses at the boundaries Sgxx = ^^^yy and 
Sayy = {2S^ + SX)eyy. Wc use a dipole approximation for 
the constraint of fixed length of the rods as earlier. The 
rigid rod is then mechanically equivalent to a force dipole 
of strength p = where e — eyy cos^ 9 and x depends on 
the linear response of the material. For rod orientations 
in a given angular range d9, the stress arising from these 
induced dipoles is given by daij = ^xea^ajii^. Thus, 



SX = !^ J cosi9fsm{9fd9 



(2S^i + SX) = 



2tt 



cos{9)'^d9 



(25) 



Solving for Sfi and 6X, we find 6^ — 5X — |x«. 

We calculate x ^ in Eq. (|4]). Here, the response 
function for an isotropic, homogeneous and compress- 
ible two dimensional system can also be written in terms 
of parallel and perpendicular components as aij{r) = 
a\\rifi + ce±iSij — fifj), with the components given by 
El: 



1 



-(1 -I- /3)(ln[r/A] + ^ie) + ln[2//3^] T ^^y^ 



(26) 



where is the Euler-Mascheroni constant. The above 
equations describe the elastic response of the material for 
distances r smaller than A, which represents an upper 
cut-off length, below which the elasticity of the material 
can be considered two-dimensional. Solving Eq. we 
find 



X = M^aV(l-f/?). 



(27) 



Together with the expressions for 5X and 5^ obtained 
above, we find the following differential equation for j3: 



dn 



(1-3/3)/? 
i(l + /3) ■ 




(28) 



FIG. 8: The flow diagram for the degree of compressibility /3 
and Poisson ratio v showing a stable fixed point a.t [3 = v = 
1/3, and an unstable fixed point at /3 = 0, = 1. 



The Poisson's ratio v, which in two dimensions is defined 
asj^ = A/(A-|-2yu), and (3 are related by the simple expres- 
sion u — \ — 2/3. Using the same framework as earlier, we 
find that at high rod densities the composite would now 
approach a steady state with a compressibility (3 — l/i 
(as earlier), but Poisson's ratio v — 1/3 as seen in Fig. [S] 
The differential equations for the Lame coefficients A and 
/i are given by : 

d/i dX ^TTo? A + 2/i 
dn dn 8 A -I- 3/i 

The elastic response in a two dimensional material has 
an upper cut-off distance, described in our case by A. 
For a 2D monolayer lying on a 3D viscous fiuid phase, 
this lengthscale is set by the ratio of the two dimensional 
shear modulus of the monolayer to the three dimensional 
modulus of the fluid phase. In this case A represents a 
crossover length below which strains are two-dimensional, 
and above which they are dominated by viscous damping 
in the three-dimensional fluid underneath the monolayer. 



VI. DISCUSSION 

We have studied the collective mechanical response of 
composites of rods embedded in an elastic medium such 
as MTs in F-actin [2^, or carbon nanotubes [3) in 
synthetic or biological gels, using a mean fleld approach 
and a dipole approximation for constraints on the rods. 
We find that an initially incompressible material, both 
in three and two dimensions, will stay incompressible 
(/3 — 0) even on adding rods to it, however, if the medium 
is even marginally compressible to begin with, the addi- 
tion of rods will drive it to a fixed point at /? = 1/3, 
signifying compressibility. In three dimensions this cor- 
responds to a stable fixed point in the Poisson's ratio at 
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1/ = l/A (compressible) and an unstable fixed point at 
v = 1/2 (incompressible) , while in two dimensions it cor- 
responds to = 1/3 and v — \ respectively- Our results 
may help to explain recent experiments 29] that have 
reported v < 1/2 for composites of microtubules and 
F-actin networks (in 3D). We also derive an expression 
for the Lame coefficients as a function of arbitrary rod 
density. Although approximate, this calculation recovers 
the expected results for very low and high rod densities. 
Our results suggest an important role for stiff filaments 
such as MTs and stress fibers in the mechanics of the 
cell cytoskeleton — they not only enhance the stiffness of 
the CTtoskeleton [4^ and its ability to bear large forces 
[30I [s^l, but may also endow it with a compressibility 
that enables it to undergo small volume changes when 
necessary. 

The ability to tune the Poisson's ratio of a composite 
by varying rod concentration has important applications 
for engineering materials. The primary mechanism of 
failure in composite materials is through a tensile fail- 
ure caused by the reinforcing fibers getting narrower and 
pulling away from the matrix when stretched. Allowing 
the composite to have a finite compressibility will lead to 
comparatively less thinning and stretching of the fibers 
for a given load, and consequently the load required to 
cause structural failure will significantly increase. Our 
result that the addition of elastic rods or fibers leads 
to a monotonic evolution of Poisson's ratio toward the 
value 1/4 in 3D, either from above or below, is further 
consistent with recent numerical calculations for fiber- 
reinforced concrete, showing a weak increase in u with 
fiber density 45[ . Although our work was primarily mo- 
tivated by intracellular networks that are nearly incom- 
pressible, concrete represents an interesting case of a ma- 
trix with Poissons ratio of ~ 0.2. While the authors of 
Ref. [45] do not make any general predictions or state- 
ments concerning Poissons ratio, they report specific val- 
ues for certain fiber volume fractions. Interestingly, they 
find that fiber inclusions lead to an increasing Poissons 
ratio (< 1/2) for all systems studied, consistent with our 
general predictions in Fig. [5] and Fig. S) Flow diagrams 
in the area Poisson's ratio have been previously reported 



for composites made of a background isotropic matrix 
and unidirectional fibers arranged randomly or in a su- 
perlattice, for t he p lane perpendicular to the long axis of 
the fibers [13, H^- These studies are applicable to 2D 
composites of circular disc-like inclusions, and report a 
fixed point in the area Poisson's ratio v ~ l/i for random 
and Kagome arrangement of discs. 

Fiber reinforced composites, such as the system we 
study, can have both direct interactions between the in- 
clusions, as well as indirect interactions through the sur- 
rounding matrix. In this manuscript we consider the 
limit of strong direct interactions between the fibers, and 
assume the fibers to be rigidly embedded in the matrix 
and employ a self-consistent approach in calculating the 
macroscopic elasticity of the composite. Our approach is 
similar in concept to homogenization methods [23] and 
differential effective medium theories ^3 that take the 
point of view that a composite material may be con- 
structed by making infinitesimal changes in an already 
existing composite. Further investigations are needed to 
take into account the effect of tuning the strength of in- 
teractions between fibers, as well as of fiber curvature 
due to the bending elasticity of the fibers and thermal 
fluctuations. It is likely that fiber curvature will lead 
to smaller effective elastic modulus in longitudinal com- 
pression than in tension in the same direction with im- 
portant consequences for the effective medium normal 
stresses. Furthermore, in the present study the distribu- 
tion of fibers has been assumed to be homogeneous and 
isotropic. In biological composites, however, one often 
comes across highly inhomogeneous regions of fiber rein- 
forcement, such as a cytoskeletal composite that shows 
bundle formation and orientational ordering of filaments. 
These features have important implications for the non- 
linear mechanical response of the cytoskeleton and will 
be addressed in future work. 
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